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基于 反应 输 运 而 合 机 理 的 腐蚀 演化 研究 
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摘 要 :油气 管道 输 运 介 质 中 化 学 反应 、 离 子 传 质 、 金 属 溶解 之 间 存 在 着 复杂 的 交互 作用 , 极 吻 诱发 
浓度 场 -电势 场 厅 合作 用 下 的 腐蚀 问题 。 本 研究 综合 考虑 溶液 化 学 反应 、 界 面 电化 学 反应 、 离 子 扩 
散 和 电 迁 移 等 多 种 不 同时 空 尺度 的 物理 化 学 过 程 ,基于 电场 浓度 场 的 耦合 机 理 建 立 腐蚀 演化 动力 
学 模型 ,分 析 了 浓度 分 布 规律 、 腐 蚀 电 流 密 度 分 布 规律 腐蚀 形 貌 演变 规律 等 动态 腐蚀 行为 ,进而 探 
究 了 局 部 化 学 环境 以 及 电化 学 行为 的 耦合 过 程 对 腐蚀 生长 动力 学 的 作用 机 制 。 研 究 发 现 腐蚀 演变 
规律 由 电学 作用 和 化 学 作用 之 间 的 竞争 作用 决定 , 随 着 时 间 的 推移 ,电学 作用 占据 主导 地 位 ,导致 
腐蚀 电流 密度 逐渐 减 小 , 蚀 坑 形 貌 的 演变 始终 保持 对 称 性 ,而 腐蚀 坑 宽 度 的 增长 小 于 深度 的 增长 。 
“关键 词 : 多 场 耦合 ;腐蚀 ;动力 学 模型 ;离子 传 质 ;演化 规律 
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Corrosion evolution based on reaction transport coupling mechanism 
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Abstract : In the transportation medium of oil-gas pipelines, there exist complex interactions among chemi- 
cal reactions , ion mass transfer and metal dissolution , which can easily induce corrosion under the coupling 
effect of concentration field and electric potential field. This paper comprehensively incorporates various 
physical and chemical processes at different time and space scales , including solution chemical reactions, 
interface electrochemical reactions , ion diffusion and electromigration. Based on the coupling mechanism of 
electric field and concentration field ,a corrosion evolution kinetic model is established. The dynamic corro- 
sion behaviors including concentration distribution, corrosion current density and corrosion morphology are 
analyzed. The mechanism of the coupling process of local chemical environment and electrochemical be- 
havior on corrosion growth kinetics is explored. This study found that the law of corrosion evolution is de- 
termined by the competition between electrical and chemical effects. As time proceeds, electrical effects 
dominate the corrosion process ,resulting in a gradual decrease in corrosion current density. The evolution 


of pit morphology always maintains symmetry. The increase in the width of the corrosion pit is less than the 
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increase in the depth. 
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随 着 C0; 驱 油 技术 的 广泛 应 用 ,油气 管道 中 的 


理 模 型 ,并 基于 腐蚀 速率 对 腐蚀 形 貌 及 计算 区 域 进 


C0, 愈 发 常见 。C0, 溶 于 水 发 生 一 系列 化 学 反应 会 
生成 H,C0; 等 物质 ,由 于 HC0; 具 有 “储存 氨 离 子 ” 
的 特性 ,导致 C0, 水 溶液 比 同样 pH 的 强酸 溶液 腐蚀 
性 更 强 趾 ,因此 CO0, 腐 蚀 问题 成 为 管道 安全 的 主要 
威胁 。 通 常 腐蚀 产物 膜 会 覆盖 在 金属 表面 ,但 是 
一 旦 产物 膜 出 现 破损 ,金属 就 会 裸露 出 来 成 为 阳极 ， 
周围 的 产物 膜 由 于 高 电位 成 为 阴极 ,最 终 形成 了 电 
化 学 反应 体系 ,在 浓度 梯度 和 电位 梯度 的 耦合 作用 
下 ,腐蚀 坑 表面 不 断 演变 ,进而 导致 管道 穿孔 等 事 
故 。 这 一 问题 已 成 为 制约 管道 安全 连续 运行 的 重大 
瓶颈 ,引起 工程 界 和 学 术 界 的 重点 关注 。 因 此 ,建立 
和 坟 耦合 作用 下 腐蚀 演化 动力 学 模型 ,探究 典型 腐 
示 况 下 的 腐蚀 演化 机 制 具 有 重要 意义 。 

轨 腐 亿 是 包含 一 系列 互相 耦合 的 物理 过 程 和 电化 
学 过 程 的 复杂 现象 ,如 传 质 过 程 .化 学 反应 .金属 / 溶 
渡 锡 面 的 电化 学 反应 。 腐 蚀 的 微小 尺度 以 及 金属 表 
面 \ 动 态 电位 浓度 梯度 之 间 的 复杂 耦合 作用 加 剧 了 
辜 记 研究 的 复杂 性 。 学 者 们 开发 了 大 量 腐蚀 模型 来 
研 索 多 因素 作用 下 的 腐蚀 演变 过 程 ,这 些 模型 可 以 
分 弦 两 类 :随机 模型 和 机 理性 模型 。 随 机 模型 主要 
找 述 大 量 腐蚀 坑 的 产生 和 生长 过 程 , 主要 包括 蒙特 
Tarea A a LIE SEU DURCH ARE 
链 人 ”等 方法 。 这 些 随机 模型 主要 描述 了 腐蚀 现象 
的 随机 性 ,并 没有 从 真实 的 物理 化 学 参数 以 及 溶液 
环 镜 因素 的 角度 去 揭示 腐蚀 机 理 。 而 机 理性 模型 主 
要 与 大 量 的 物理 化 学 参数 紧密 联系 在 一 起 ,采用 多 
物理 场 厢 全 方法 ,基于 物质 传输 定律 和 电化 学 动力 
学 定律 ,用 物理 模型 和 数学 方程 来 传递 出 腐蚀 坑 的 
形 貌 演变 规律 以 及 腐蚀 系统 内 的 浓度 场 .电势 场 分 
布 规律 。 学 者 们 提出 了 很 多 机 理性 模型 ”来 描述 
腐蚀 坑内 的 传 质 过 程 和 化 学 反应 ,并 且 考 虑 了 2 种 
反应 对 传 质 过 程 的 影响 :金属 /溶液 界面 的 电化 学 
反应 ;@ 腐 蚀 溶 液 中 的 化 学 反应 。 腐 蚀 形 貌 随时 间 
的 演变 规律 是 影响 腐 亿 过 程 的 关键 因素 ,而 大 多 数 
腐蚀 模型 并 未 考虑 腐蚀 形 貌 的 动态 演变 过 程 对 腐蚀 
行为 的 影响 。 此 外 ,目前 的 腐蚀 模型 中 很 少 考虑 电 
极 反应 、 化 学 反应 与 多 组 分 传 质 过 程 之 间 的 耦合 作 
用 。 因 此 ,仍然 需要 基于 化 学 .电化 学 、 动 电学 和 传 
质 学 的 耦合 机 理发 展 一 个 综合 性 的 多 物理 场 耦合 机 


n 


行 实时 动态 更 新 ,捕捉 腐蚀 的 动态 演变 过 程 ,以 反映 
更 真实 的 腐蚀 演化 过 程 。 

本 研究 建立 浓度 场 -电场 耦合 作用 下 的 腐蚀 动 
力学 模型 ,分 析 了 浓度 分 布 规律 ,腐蚀 电流 密度 分 布 
规律 腐蚀 形 貌 演变 规律 等 动态 腐 介 行为 ,进而 探究 
局 部 化 学 环境 以 及 电化 学 行为 的 耦合 过 程 对 腐蚀 生 
长 动力 学 的 作用 机 制 ,并 建立 腐蚀 行为 与 腐蚀 坑内 
部 电化 学 环境 的 内 在 联系 ,从 化 学 和 电学 因素 耦合 
作用 的 角度 揭示 腐蚀 演化 规律 。 


1 CO; 腐蚀 体系 


本 研究 以 近 管 壁 区 包含 C0, 与 游离 水 的 酸性 溶 
液体 系 为 研究 对 象 , 主要 考虑 C0, 溶 于 水 的 四 步 平 
衡 反 应 以 及 Fe ' 水解 反 应 ,建立 基于 离子 反应 输 运 
机 理 的 金属 管 壁 腐蚀 机 理 模 型 。 通 常 油气 管线 中 
CO, 的 分 压 不 可 忽略 ,而 且 近 管 壁 区 溶液 的 pH 值 也 
相对 较 小 ,此 时 可 忽略 腐蚀 电化 学 反应 中 H,O 的 还 
原 反 应 ,而 考虑 Fe 的 阳极 反应 以 及 H* 、H,C03 还 原 
的 阴极 反应 。 因 此 ,整个 反应 体系 中 化 学 反应 和 电 
化 学 反应 如 下 。 


1 ) 化 学 反应 过 程 
CO Fk 
CO, (g)===CO, (aq) (1) 
CO, 水 合作 用 
C0,(g) + H,O ——CO, (aq) (2) 
H,C0; 部 分 电离 
H,CO, —H* «HCO; (3) 
HCO; 部 分 电离 
HCO; —5H* «CO?- (4) 
Fe “水解 反 应 
Fe'* + H,O 一 一 Fe(OH) * -H* (5) 
2) 电 化 学 反应 过 程 
阳极 反应 
Fe —Fe* +2e- (6) 
阴极 反应 主要 是 也 ”的 还 原 过 程 
2H* +2e —9H, (7) 


在 管 壁 附近 还 会 发 生男 一 个 重要 的 反应 ， 
H,C0; 的 直接 还 原 反 应 , 即 
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2H,CO, 42e; —9H, +2HCO7 (8) 

腐蚀 模型 示意 图 如 图 1 所 示 , 其 中 上 边界 为 开 
边界 (主体 溶液 ) ,下 边界 为 管 壁 ,左右 边界 为 电解 
液 边界 ,假设 初始 的 腐蚀 坑 为 半圆 形 。Mohammed 
A 7 发 现 管线 钢 在 C0, 溶 液 中 的 点 蚀 坑 尺寸 分 布 
范围 在 9 ~ 100 pm 之 间 , 因 此 本 研究 选择 处 于 该 范 
围 中 的 50 pm 作为 点 蚀 坑 初 始 半径 进行 数值 模拟 。 


CO, 主体 溶液 


i 
CO,*H,O— H.CO, HCO; +H 
i 


扩散 & 电 迁移 C0 


H,CO,*e “H+HCO, Fe—Fe”+2e H+e=H 
x |(0, 0) EY Fe'-H,O—FeOH *H 


阴极 表面 阴极 表面 
阳极 点 蚀 坑 


图 1 点 蚀 系统 模型 示意 图 


Fig.1 Schematic diagram of pitting system model 


多 场 耦合 腐蚀 动力 学 模型 


二 本 研究 所 建立 的 腐蚀 演化 动力 学 模型 主要 基于 
浓 麻 场 . 电 势 场 的 多 物理 场 而 合 特性 ,并 考虑 了 腐蚀 
形 训 的 动态 追踪 ,由 3 部 分 组 成 : 传 质 模型 .电学 模 
TE RUE AGE E UI 


传 质 模 型 


局 该 模型 综合 考虑 了 浓度 梯度 驱动 的 扩散 传 质 和 
由 昌 场 导致 的 电 迁 移 传 质 。 
斧 根 据 不 同 离子 的 质量 守恒 
Q c; 

m eM -R, (9) 
其 中 :c， 和 N; 分 别 为 第 i 种 离子 的 浓度 与 通 量 ;R， 
为 物质 i 的 反应 源 项 , 且 
F Reo, 


r-1 0 0 0 
[o 1 -1 0 Okceo,— henco, 
Ruco; 0 1 -1 0 k, caCHco, = k, caCu* Cuco; 
Roo: 0 0 1 0|‖ 不 acaco: — Es Cu Ccor- 
Ron 0 0 0 l kj ua = k, Cu: Con- 
Ry 0 : Log keeps = hs jch Cre con): 
L Ryo - 0 0 0 0 
(10) 
式 中 的 等 反应 常数 可 查阅 文献 [21 ] 获得。 
由 能 斯 特 - 普 朗 克 方 程 可 得 
N, = - D, Vc; - au, Fc; Vo (11) 


ES, E: ETUA A a Ch RC LE Cl 


hinay i&d HHT 
IInaA IVE 1FEH 


889 
其 中 :等 号 右边 第 一 项 为 扩散 传 质 作用 ,第 二 项 为 电 
迁移 传 质 作用 ;D; 为 组 分 i 的 扩散 系数 ;zi 是 组 分 i 的 
电荷 数 ;F 是 法 拉 第 常数 (96 500 C mol) ;9 是 电 
位 ,可 由 后 面 的 动 电学 模型 计算 得 出 ;uw 是 组 分 i 
的 物质 迁移 率 , 且 


D, 
"mi RT 
其 中 :RR 是 通用 气体 常数 (8.314/J * mol - K);T 
是 绝对 温度 。 
对 于 不 参加 电极 反应 的 物质 
ðc; 


u 


(12) 


对 于 参加 电极 反应 的 物质 
-p Lo. (14) 


2.2 电学 模型 


计算 区 域内 电位 分 布 须 满足 以 下 泊 松 方程 , 即 


Ve = -一 (15) 
其 中 :9 为 电位 ;F 为 法 拉 第 常数 ;se 为 电解 液 介 电 
常数 。 
2.3 ”腐蚀 形 貌 追踪 模型 
本 研究 将 拉 普 拉 斯 平滑 法 嵌 套 到 数值 模型 中 作 
为 任意 拉 格 朗 日 - 欧 拉 方 法 的 网 格 更 新 方法 来 保证 
网 格 质 量 , 可 用 数学 方程 表示 如 下 。 


0 x 0 yx 
REESE «S =0 16 
aX ot Y ðt (16) 
2 2 P 
P y IY y LE ap (17) 


aX! ðt oY ðt 
其 中 :x 和 y 是 空间 框架 下 的 空间 坐标 ;X 和 了 是 材 
料 框 架 下 的 参考 坐标 。 
壁面 移动 速度 是 由 金属 溶解 造成 的 , 可 由 法 拉 
第 定律 计算 , 即 


EN as 
"ongPpyg. 

其 中 :凡是 壁面 移动 的 法 向 速度 ;mr* 是 阳极 反应 的 
电荷 转移 数 ;Wr 和 pr 分 别 是 溶解 金属 Fe 的 摩尔 质 
量 和 密度 。 

在 任意 拉 格 朗 日 - 欧 拉 方 法 中 ,如果 网 格 变形 过 
大 会 导致 网 格 质量 下 降 , 计算 过 程 中 可 能 会 频繁 出 
错 ,为 了 解决 这 一 问题 ,采用 将 计算 区 域 重 新 剖 分 ， 
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在 网 格 质量 严重 退化 之 前 将 变量 的 解 从 旧 网 格 映射 
到 新 网 格 。 

2.4 ”腐蚀 产物 膜 影响 下 的 传 质 修正 


在 C0, 腐 蚀 过 程 中 , 当 Fe 和 COS. 超出 溶液 
中 的 溶解 度 极限 后 就 会 沉淀 下 来 生成 FeCO, , 即 
Fe'* + CO]. =—FeCO, (s) (19) 
24 FeCO, DCUE TE 5] e E HT , S t — BB 2) BT 
表面 从 而 阻止 离子 参与 腐蚀 过 程 ,因此 Feco, f C 
淀 有 可 能 会 阻碍 腐蚀 过 程 。 
沉积 速度 可 利用 下 式 '“! 计 算 。 


E 
R recos 三 e” R y EG -1) (20) 
FeC0; 的 过 饱和 度 S 定义 为 
> s Cre2+ co? - 21 
( TUE (21) 


DOFeC0, 的 溶解 度 极限 K TI Hr FR 49 ERR 


d —59.349 8 0.041377 x T - ^7 +24. 572 4 x logT «2.518 x ^ -0.657 x I 
= 10 


(22) 
HE 是 离子 强度 。 
qui Renan 质 特 性 对 扩散 系数 作 如 
RBE”, 


D, 
D, ;7D(l-ep) tq —— Ef (23) 
A, PE 


jb. D, ,为 组 分 ; 在 产物 膜 中 的 扩散 系数 ; D; 为 组 
Ari 在 主流 溶液 中 的 扩散 系数 ;Nu pe 为 MaeMullin 数 
VES 


Nume = (24) 
这 
其 中 :7 是 多 孔 介 质 的 弯曲 度 ;ej 为 多 孔 介 质 的 孔 院 
率 ,分 别 可 表示 为 
peu (25) 


e, =580 xe 9997-1302? £0.06 — (26) 


同 理 , 物 质 的 迁移 率 也 需要 进行 相应 的 修正 。 
D. 
wy = RY (27) 
由 于 本 研究 工 况 中 FeC0; 的 过 饱和 度 较 低 , 因 
此 针对 生成 产物 膜 的 工 况 ,只 考虑 其 产物 膜 对 传 质 
的 影响 ,忽略 其 他 影响 。 
以 上 腐蚀 动力 学 模型 的 边界 条 件 如 表 1 所 示 ， 
Con ii 是 阳极 反应 电 


学 报 b hi n ax v TER HJ H JF I) 第 39 卷 
表 1 腐蚀 动力 学 模型 边界 条 件 


Tab.1 Boundary conditions of corrosion dynamics model 


A - 形 貌 追踪 
边界 传 质 模型 电学 模型 
模型 
上 边界 ci 6j, =g 99 29 dx 2dyz0 
on 
左右 
-n:N;z0 sg dep dx =dy=0 
边界 on 
阳极 N; Viyipe g 99. E Ys V, - ly My, 
边界 nF on 7 npe FPre 
阴极 
ne N; = Ë 2g99 à Xs dy 20 
力 界 nF on 


在 模拟 计算 过 程 中 ,首先 根据 动 电 学 模型 计算 
腐蚀 系统 的 电位 分 布 特性 ,然后 根据 传 质 模型 计算 
瞬 态 浓度 场 分 布 规 律 。 考 虑 到 计算 稳定 性 和 收敛 性 
的 需求 ,本 研究 的 多 物理 场 问题 的 求解 采用 全 耦合 
方式 :所 有 的 单元 方程 首先 被 装配 成 全 局 方程 ,然后 
获得 一 个 包含 所 有 求解 变量 (电位 和 浓度 ) 的 全 局 
和 矩阵 ,对 该 全 局 矩阵 进行 求解 可 同时 获得 所 有 变量 
的 解 。 在 此 基础 上 ,基于 法 拉 第 定律 计算 腐蚀 生长 
速率 。 同 时 ,采用 任意 拉 格 朗 日 - 欧 拉 方法 ,基于 腐 
刨 生长 速度 来 追踪 腐蚀 形 貌 的 演变 过 程 并 更 新 计算 
区 域 。 最 终 可 以 捕捉 到 任意 时 刻 和 任意 位 置 的 浓度 
场 . 电 势 场 和 腐蚀 形 貌 等 腐蚀 特征 。 


3 ”模拟 结果 及 讨论 


3.1 模型 验证 


为 了 验证 本 研究 所 建立 的 数学 模型 的 准确 性 
模拟 了 pH 6 环境 中 的 腐蚀 电流 随时 间 的 演变 情况 ， 
然后 将 Zhang 等 的 同一 工 况 下 的 大 量 实验 数据 
进行 平均 化 处 理 , 再 将 其 实验 结果 与 本 研究 的 模型 
预测 结果 进行 对 比 。 如 图 2 所 示 , 随 着 时 间 的 推移 ， 
预测 腐蚀 电流 与 实验 腐蚀 电流 均 逐 渐 减 小 直至 稳 
定 , 这 是 由 于 腐蚀 产物 膜 逐 渐 生 长 阻碍 了 腐蚀 过 程 ， 
由 图 2 可 知 本 研究 的 模拟 结果 与 实验 结果 吻合 较 
好 ,并 且 准 确 地 反映 了 时 间 对 腐蚀 行为 的 影响 ,证 明 
了 本 研究 所 建立 的 电场 -浓度 场 耦合 作用 下 CO, JA 
蚀 动 力学 模型 的 准确 性 。 
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图 2 腐蚀 电流 演变 过 程 的 模拟 与 实验 对 比 


Fig.2 Comparison of the experimental and calculated 


galvanic current evolution 


3.2 浓度 场 特征 分 析 


腐蚀 体系 内 不 同 物质 的 浓度 分 布 是 传 质 

化 学 反应 .电极 反应 综合 作用 的 dcs 
过 各 的 耦合 作用 。 图 3 为 由 管 a 
Cafe tte ) D s a T KRE 58 EAW YR RE 2 
IÈI c; - co) 分 布 规律 。 由 图 3 T Fe^* 浓度 在 
区 名 坑内 部 最 高 ,这 是 由 于 金属 溶解 发 生 在 腐蚀 坑 
HG JEPEN Pe ,然后 Fe 自 下 向 上 从 腐 人 
境 禾 口 流出 进入 主流 溶液 ;同时 H, C0; 还 原生 成 
HEQ; ,由 于 H,C0, 在 金属 表面 被 消耗 ,C0, 需 不 断 
地 向 金属 表面 运 移 以 发 生 水 合 反 应 生成 H CO, P 
HEGO, KYK BE rh E BEYA EEE T E 系 渐 减 小 ,从 而 反 
喘 俐 其 传输 方向 为 从 主流 到 管 壁 。C0, 和 HCO; 均 
珊 出 与 主流 溶液 浓度 的 较 大 浓度 偏差， 显示 了 
CF 冰 合 反应 的 显著 影响 。 腐 蚀 所 需要 的 HUI 
H4€0, ,主要 是 主流 向 壁面 的 传 质 过 程 所 提供 的 ,由 
alt E i KI C0, 的 水 合 反应 速率 很 慢 ,不 足以 及 
时 观 供 腐蚀 所 需要 的 H* 和 HS CO. ,但 可 作为 腐蚀 过 
程 中 H* 和 HC0; 的 重要 补充 来 源 。 值 得 注意 的 是 ， 
由 于 本 研究 所 模拟 的 溶液 为 近 中 性 环境 ,HH’ 浓度 的 
变化 量 级 约 为 10 mol * m^? «10 ? mol - m ,而 其 
他 离子 的 浓度 变化 量 级 约 为 10-!'mol * m ,因此 虽 
然 实 际 上 图 3 中 也 + 浓度 发 生 了 变化 ,但 是 看 起 来 并 
不 明显 ,可 结合 图 7 看 出 H' 浓度 在 不 同 区 域 的 显著 
变化 。 


3.3 ”腐蚀 形 貌 演变 规律 


图 4 所 示 为 腐蚀 坑 形 貌 随时 间 的 演变 规律 , 腐 
刨 坑 形 貌 的 演变 始终 保持 对 称 性 。 此 外 ,腐蚀 坑 的 
深度 和 宽度 均 随 着 时 间 推 移 而 逐渐 增长 ,40h 后 , 腐 
蚀 坑 的 深度 由 最 初 的 50. um 演化 为 53. 3 pum, 宽度 
由 最 初 的 100 um 演化 为 105.5 pm。 图 5 定量 比较 


基于 反应 输 运 耦合 机 理 的 腐蚀 演化 研究 ~ 
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了 腐蚀 坑 深度 和 宽度 随时 间 的 变化 规律 ,可 以 发 现 
窝 刨 坑 宽度 的 增长 略 小 于 深度 的 增长 ,这 一 现象 可 
以 通过 后 文 分 析 腐 蚀 坑内 腐蚀 电流 密度 的 分 布 规律 


c;-c, (mol * m?) 


02L LLLE eo e © 
-60 -40 20 0 20 40 60 480 100 
y/um 


图 3 垂直 管 壁 方向 各 物质 浓度 与 主流 
溶液 浓度 的 差 值 分 布 规律 
Fig.3 Deviation of species concentrations from the bulk values 


along the vertical direction away from the metal surface 


200 220 240 260 280 300 320 
x/um 


图 4 不 同时 刻 腐蚀 坑 形 貌 演变 规律 


Fig.4 Time-dependent deformation of pit geometry 


Wh 


图 5 不 同时 刻 腐蚀 坑 宽 度 与 深度 的 对 比 
Fig.5 Time-dependent change of pit depth and width 


3.4 腐蚀 电流 密度 变化 规律 


如 图 6 所 示 ,腐蚀 坑 侧 壁 的 腐蚀 电流 密度 明显 
小 于 底部 的 腐蚀 电流 密度 ,因此 随 着 时 间 的 推移 , 腐 
蚀 坑 宽 度 的 增长 会 逐渐 小 于 深度 的 增长 ,同样 导致 
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腐蚀 朝 着 管 壁 更 深 处 发 展 , 从 而 引发 管道 穿孔 。 值 
得 注意 的 是 ,在 腐蚀 演化 过 程 中 ,腐蚀 电流 密度 随 着 
时 间 推 移 逐 渐 减 小 ,这 一 现象 同样 可 以 通过 分 析 化 
学 作用 与 电学 作用 的 耦合 效应 来 解释 。 


x/um 


图 6 不 同时 刻 腐蚀 坑 腐 蚀 电流 密度 分 布 规律 

T=  Fig.6 Distribution of corrosion current density 

> along the pit wall at different time 
钙化 学 作用 如 图 7 和 图 8 所 示 。 图 7 显示 了 不 同 
时 到 H 浓度 沿 着 腐 他 坑 壁 面 的 分 布 规律 ,可 以 发 现 
$ 鲁 坑 两 侧 靠 近 阴 极 区 域 的 H 浓度 明显 较 低 。 这 
是 因为 阴极 区 域 的 还 原 反应 消耗 了 大 量 H* ,造成 明 
极 表 面 与 阳极 点 蚀 坑 交界 处 的 浓度 梯度 较 大 ,促进 
抠 间 蚀 坑 两 侧 靠近 阴极 区 域 的 H+* 向 附近 阴极 区 域 
护 散 ,同时 阴极 表面 与 阳极 点 刨 坑 交界 处 的 电位 差 
世 潮 成 了 点 他 坑 两 侧 FL 向 附近 阴极 表面 的 电 迁 移 
传 质 。 因 此 ,点 包 坑 两 侧 靠近 阴极 区 域 的 H' 浓度 较 
信 忆 同时 ,H* 浓度 随 着 时 间 的 增加 而 逐渐 增 大 。 交 
ML UR HE 加 是 H' 浓度 的 函数 ,可 以 反映 出 H' 浓 
度 对 腐蚀 电流 密度 的 影响 ,因此 图 8 显示 交换 电流 
ER 
Af FEL; 5 EONI E H * 浓度 随时 间 的 增 
长 趋势 造成 的 。 
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图 7 不 同时 刻 腐蚀 坑内 HL "浓度 分 布 规律 


Fig.7 Distribution of H* concentration along 


300 310 


the pit wall at different time 
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i/(107A * m?) 


器 


8 不 同时 刻 腐 蚀 坑内 交换 电流 密度 分 布 规律 


Fig.8 Distribution of exchange current density 


along the pit wall at different time 

电学 作用 如 图 9 和 图 10 所 示 。 图 9 显示 了 不 
同时 刻 沿 着 腐蚀 坑 壁 面 的 过 电位 分 布 规律 ,首先 可 
以 发 现 过 电位 的 分 布 同样 呈现 对 称 性 的 特征 ,此 外 ， 
过 电位 随 着 时 间 的 推移 逐渐 减 小 , 即 电学 效应 对 腐 
蚀 演 化 的 驱动 作用 逐渐 减弱 ,因此 图 10 显示 过 电位 
对 腐蚀 电流 密度 的 影响 i, 随 时 间 推 移 逐 渐 减 小 , 且 
同样 呈现 对 称 分 布 的 特点 。 


77.4 10h 


x/um 


图 9 不 同时 刻 腐 刨 坑 表 面 过 电位 分 布 规律 
Fig.9 Distribution ofoverpotential along 


the pit wall at different time 


86 x 10h » 
85 
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x/um 
图 10 不 同时 刻 腐蚀 坑 表 面 i, 分 布 规律 


Fig. 10 Distribution of i, along the pit wall at different time 


综 上 所 述 , 随 着 时 间 的 推移 ,溶液 的 侵蚀 性 越 来 
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越 强 ,但 是 过 电位 迅速 减 小 ,导致 i, 的 逐渐 增 大 以 及 
,的 逐渐 减 小 。 这 种 化 学 作用 和 电学 作用 之 间 的 竞 
争 效应 决定 了 腐蚀 电流 密度 的 分 布 规律 及 发 展 趋 
势 。 由 于 腐蚀 电流 密度 是 i 和 i, 的 乘积 ,对 称 的 io 
分 布 轮廓 和 对 称 的 i, 分 布 轮廓 共同 导致 了 腐蚀 电流 
密度 分 布 的 对 称 性 。 随 着 时 间 的 推移 ,i, 增 大 的 程 
度 明 显 小 于 i, 减 小 的 程度 , 即 电学 作用 比 化 学 作用 
更 重要 ,导致 腐蚀 电流 密度 逐渐 减 小 。 

此 外 ,为 了 对 比分 析 腐蚀 产物 膜 对 于 腐蚀 行 > 
的 影响 ,图 11 显示 了 腐蚀 产物 膜 存在 与 去 除 腐蚀 产 
物 膜 条 件 下 的 点 蚀 坑 表 面 腐蚀 电流 密度 分 布 规律 。 
可 以 发 现 腐蚀 产物 膜 存在 时 的 腐蚀 电流 密度 明显 小 
于 无 腐蚀 产物 膜 条 件 下 的 腐蚀 电流 密度 ,这 是 由 于 
腐蚀 产物 膜 内 部 的 多 孔 介质 结构 阻碍 了 多 种 组 分 的 
传 质 过 程 ,从 而 减缓 了 一 系列 化 学 反应 和 电化 学 反 
庶 是 必需 的 反应 物 /生成 物 运 移行 为 ,最 终 降 低 了 腐 
蚀 电 化 学 反应 速率 。 


~ 


Q2 


腐蚀 产物 膜 存在 /去 除 时 的 点 蚀 坑 腐 蚀 电流 密度 对 比 
Fig. 11 


the presence/removal of corrosion product film 


Corrosion current density comparison with 


4 结论 


本 研究 针对 化 学 反应 、 离 子 传 质 、 电 极 界面 反应 
综合 作用 下 的 腐蚀 问题 ,基于 浓度 场 与 电势 场 的 耦 
A BUE ,建立 腐蚀 演化 动力 学 模型 ,探究 了 浓度 场 、 
腐蚀 形 狐 .腐蚀 电流 密度 随时 间 的 变化 规律 ,重点 从 
化 学 和 电学 因素 耦合 作用 的 角度 揭示 腐蚀 演化 规 
律 。 主 要 结论 如 下 。 

1) 本 研究 建立 的 多 场 耦合 模型 所 预测 的 腐蚀 电 
流 与 实验 平均 结果 的 吻合 度 较 好 ,并 能 准确 反映 出 
腐蚀 演化 过 程 的 时 间 效 应 ,证 实 了 腐蚀 演化 动力 学 
模型 的 可 靠 性 。 

2) Fe * HCO; 作为 反应 产物 ,其 浓度 由 金属 辟 
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面 到 主体 溶液 逐渐 降低 ,而 C0, 由 主体 溶液 向 壁面 
运 移 以 发 生 水 合 反 应 生成 不 断 被 消耗 的 H,C0,。 

3) 由 于 浓度 分 布 以 及 电位 分 布 的 对 称 性 ,腐蚀 
坑 形 貌 的 演变 始终 保持 对 称 性 ,而 腐蚀 坑 宽 度 的 增 
长 略 小 于 深度 的 增长 。 

4) 对 称 的 记分 布 轮 廊 和 对 称 的 i, 分 布 轮廓 共同 
导致 了 腐蚀 电流 密度 分 布 的 对 称 性 ; 随 着 时 间 的 推 
移 , 电 学 作用 占据 主导 , 导致 腐蚀 电流 密度 逐渐 
减 小 。 


参考 文献 : 


[1] NESIC S,LI H, HUANG J,et al. An open source mechanistic model 
for CO,/H,S corrosion of carbon steel[ C ]//NACE International. 


Proceeding of 2009 Corrosion Conference and Expo. Atlanta; 
OnePetro ,2009 :09572. 

[2] 李 佳 航 , 王 丹 , 谢 飞 ,等 . 油气 管道 的 COS 腐蚀 及 防护 研究 进展 
[J]. 表面 技术 ,2021,50(4) :172-183. 
LI Jiahang, WANG Dan, XIE Fei,et al. Research progress on CO, 


corrosion and protection of oil and gas pipelines[ J |. Surface tech- 
nology ,2021 ,50(4) :172-183 (in Chinese). 

[3] YUE X Q, REN Y Q, HUANG L Y et al. The role of C17 in the 
formation of the corrosion products and localised corrosion of 15Cr 
martensite stainless steel under an CO; -containing extreme oilfield 
condition[ J ] . Corrosion science ,2022 ,194 109935. 

[4] 骆 正 山 , 赵 乐 新 , 王 小 完 . 基于 动态 贝 叶 斯 网 络 的 海底 管道 点 
蚀 疲 劳损 伤 失 效 模 型 研究 [J]. 表面 技术 ,2020,49 (1): 
269-275. 

LUO Zhengshan,ZHAO Lexin, WANG Xiaowan. Failure model for 


pitting fatigue damaged pipeline of subsea based on dynamic bayes- 
ian network [ J ]. Surface technology, 2020,49 (1) : 269275 ( in 
Chinese). 

[5] OSSAIC I, BOSWELL B,DAVIES I J. Predictive modelling of in- 
ternal pitting corrosion of aged non-piggable pipelines[ J]. Journal 
of the electrochemical society ,2015 ,162 (6) :C251-C259. 

[6] GUO Z W,CHEN H. Monte carlo simulation method for pit erosion 
process of cable wire based on bi-plane cellular automata| J ] . Inter- 
national journal of steel structures ,2020 ,20 :1527-1539. 

[7] TRH, EFA, EHN, S. 基于 元 胞 自动 机 的 飞机 蒙 皮 表 面 
点 刨 过 程 模拟 [ 1] - 腐蚀 与 防护 ,2019 ,40(6) :55-61. 
DING Qingmiao, WANG Yujun, CUI Yanyu, et al. Simulation of 


pitting corrosion of aircraft surface based on cellular automat [ J ]. 
Corrosion & protection ,2019 ,40 (6) :55-61 (in Chinese). 

[8] 王 宁 .基于 元 胞 自动 机 法 的 金属 局 部 腐蚀 损伤 演化 研究 [D1]. 
天 津 :中 国民 航 大 学 ,2020. 

[9] XUZH,LU J F, WEI X L,et al. 2D and 3D cellular automata sim- 


ulation on the corrosion behaviour of Ni-based alloy in ternary mol- 
ten salt of NaCI—KCI—ZnCl, [ J]. Solar energy materials and solar 
cells ,2022 ,240 111694. 

[10] KER, BES E, A. 元 胞 自动 机 在 金属 材料 腐蚀 研究 中 


微 信 公众 号 :应 用 力学 学 报 


894 


的 应 用 [中 .材料 导报 ,2023 ,37 (8) :21050078. 
ZHANG Xiqing, TENG Yingxue, GUO Jing, et al. Application of 


cellular automata in the research of metal material corrosion [ J ]. 
Materials reports ,2023 ,37(8) :21050078 (in Chinese). 
[11] OSSAI C I, BOSWELL B, DAVIES I. Markov chain modelling for 
time evolution of internal pitting corrosion distribution of oil and gas 
pipelines[ J ] . Engineering failure analysis ,2016 ,60 209-228. 
[12] WANG Y,ZHANG P,QIN G. Reliability assessment of pitting cor- 
rosion of pipeline under spatiotemporal earthquake including spa- 
tial-dependent corrosion growth[ J ] . Process safety and environmen- 
tal protection ,2021 ,148 1166-178. 
[13] LIU C, KELLY R G. A review of the application of finite element 
method ( FEM) to localized corrosion modeling [ J ]. Corrosion, 
2019 ,75:1285-1299. 
[14] JAFARZADEH S, CHEN Z, BOBARU F. Peridynamic modeling of 
repassivation in pitting corrosion of stainless steel[ J]. Corrosion, 


372018 ,74 :393-414. 
[ 15> WANG K,LI C,LI Y,et al. A fully coupled model of hydrodynam- 


Oic chemical electrochemical processes for CO, uniform corrosion in 
multi-physiecs environment [ J ]. Journal of petroleum science and 
© insi ,2020,193 107436. 
[ e SAEEDIKHANI M, WIJESINGHE S, BLACKWOOD D J, Moving 
"boundary simulation and mechanistic studies of the electrochemical 
X md protection by a damaged zinc coating[ J ]. Corrosion sci- 
Cf ence ,2020,163 : 108296. 
[EN ONISHI Y,TAKIYASU J, AMAYA K,et al. Numerical method for 
O time-dependent localized corrosion analysis with moving boundaries 
5 .by combining the finite volume method and voxel method| J]. Cor- 
5 rosion science ,2012 ,63 :210-224. 
ik MAI W,SOGHRATI S. New phase field model for simulating gal- 
vanic and pitting corrosion processes [ J|. Electrochimica acta, 
(—2018 ,260 :290-304. 
[ 1917 GAO H D, JU L L,LI X ,et al. A space-time adaptive finite element 
"Toa with exponential time integrator for the phase field model of 
pitting corrosion [ J ]. Journal of computational physics, 2020, 
406 :109191. 
[20] MOHAMMED S A,HUA Y,BARKER R ,et al. Effect of calcium 
on X65 carbon steel pitting in saturated CO, environment [ J ]. 
Electrochimica acta ,2022 ,407 ; 139899. 
[21] PALMER D A, VAN E R. The chemistry of metal carbonate and 
carbon dioxide complexes [ J ]. Chemical reviews, 1983,83 (6) : 
651-731. 
[22] JOHNSON M,TOMSON M. Ferrous carbonate precipitation kinetics 
and its impact CO2 corrosion [| C ]//NACE International. Proceed- 


ings of 1996 Corrosion Conference and Expo. Houston ; OnePetro, 
1996 :268. 


[23] SUN W,NESIC S, WOOLLAM R C. The effect of temperature and 


投稿 网 站 : http :// cjam. xjtu. edu. en 


应 用 力学 学 报 


[24] 


[25] 


ChinaXiv 合 作 期 刊 第 39 郑 


ionic strength on iron carbonate ( FeCO, ) solubility limit[ J ]. Cor- 
rosion science,2009,51(6) 1273-1276. 

HI C. 腐蚀 产物 沉积 对 局 部 腐蚀 行为 影响 的 数值 模拟 L[D]. 大 
连 :大 连理 工大 学 ,2013. 

ZHANG G A,CHENG Y F. Localized corrosion of carbon steel in a 


CO; -saturated oilfield formation water[ J ]. Electrochimica acta, 
2011,56 (3) :1676-1685. 
(编辑 WE ERR) 


TY 


qon dn 


BÀ RA A 


MBPS EAR, 1970 年 生 ， 山 东 济 南 
市 人 ， 杨 树 美育 法 创始 者 、 独 立 文史 艺术 学 
者 ; 季羡林 国际 文化 书院 研究 员 、 季 美 林 
国际 文化 书院 书画 委员 、 书 画家 。 


> a 
< 人 << < 


